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We show how the prevailing majority opinion in a population can be rapidly reversed by a small 
fraction p of randomly distributed committed agents who consistently proselytize the opposing opin- 
ion and are immune to influence. Specifically, we show that when the committed fraction grows 
beyond a critical value pe ~ 10%, there is a dramatic decrease in the time, Te, taken for the entire 
population to adopt the committed opinion. In particular, for complete graphs we show that when 
p < pe, Te ~ exp(a(p)N), while for p > pe, Te ~ InN. We conclude with simulation results for 
Erdos-Rényi random graphs and scale-free networks which show qualitatively similar behavior. 


PACS numbers: 87.23.Ge, 89.75.Fb 
I. INTRODUCTION 


Human behavior is profoundly affected by the influ- 
enceability of individuals and the social networks that 
link them together. Well before the proliferation of online 
social networking, offline or interpersonal social networks 
have been acknowledged as a major factor in determining 
how societies move towards consensus in the adoption of 
ideologies, traditions and attitudes ji g. As a result, 
the dynamics of social influence has been heavily studied 
in sociological, physics and computer science literature 

In the sociological context, work on diffusion of 
innovations has emphasized how individuals adopt new 
states in behavior, opinion or consumption through the 
influence of their neighbors. Commonly used models for 
this process include the threshold model [s] and the Bass 
model (9). A key feature in both these models is that once 
an individual adopts the new state, his state remains un- 
changed at all subsequent times. While appropriate for 
modeling the diffusion of innovation where investment in 
a new idea comes at a cost, these models are less suited 
to studying the dynamics of competing opinions where 
switching one’s state has little overhead. 


Here we address the latter case. From among the vast 
repertoire of models in statistical physics and mathemati- 
cal sociology, we focus on one which is a 2-opinion variant 
of the Naming Game (NG) and that we re- 
fer to as the binary agreement model. The evolution of 
the system in this model takes place through the usual 
NG dynamics, wherein at each simulation time step, a 
randomly chosen speaker voices a random opinion from 
his list to a randomly chosen neighbor, designated the 
listener. If the listener has the spoken opinion in his list, 
both speaker and listener retain only that opinion, else 
the listener adds the spoken opinion to his list (see Ta- 
ble 1). The order of selecting speakers and listeners is 
known to influence the dynamics, and we stick to choos- 
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ing the speaker first, followed by the listener. 

It serves to point out that an important difference 
between the binary agreement model and the predomi- 
nantly used opinion dynamics models Z! Bidia is that 
an agent is allowed to possess both opinions simultane- 
ously in the former, and this significantly alters the time 
required to attain consensus starting from uniform initial 
conditions. Numerical studies in have shown that for 
the binary agreement model on a complete graph, start- 
ing from an initial condition where each agent randomly 
adopts one of the two opinions with equal probability, 
the system reaches consensus in time Te ~ In N (in con- 
trast, for example, to Te ~ N for the voter model). Here 
N is the number of nodes in the network, and unit time 
consists of N speaker-listener interactions. The binary 
agreement model is well suited to understanding how 
opinions, perceptions or behaviors of individuals are al- 
tered through social interactions specifically in situations 
where the cost associated with changing one’s opinion is 
low [19], or where changes in state are not deliberate 
or calculated, but unconscious ld. Furthermore, by its 
very definition, the binary agreement model may be ap- 
plicable to situations where agents while trying to influ- 
ence others, simultaneously also have a desire to reach 
global consensus lil. 

Another merit of the binary agreement dynamics in 
modeling social opinion change seems worth mentioning. 
Two state epidemic-like models of social “contagion” (ex- 
amples in BA) suffer from the drawback that the rules 
governing the conversion of a node from a given state 
to the other are not symmetric for the two states. In 
contrast, in the binary agreement model, both singular 
opinion states are treated symmetrically in their suscep- 
tibility to change. 

Here we study the evolution of opinions in the binary 
agreement model starting from an initial state where all 
agents adopt a given opinion B, except for a finite frac- 
tion p of the total number of agents who are committed 
agents and have state A. Committed agents, introduced 
previously in R3, are defined as nodes that can influ- 
ence other nodes to alter their state through the usual 
prescribed rules, but which themselves are immune to 


influence. In the presence of committed agents adopting 
state A, the only absorbing fixed point of the system is 
the consensus state where all influenceable nodes adopt 
opinion A - the opinion of the committed agents. The 
question that we specifically ask is: how does the con- 
sensus time vary with the size of the committed fraction? 
More generally, our work addresses the conditions under 
which an inflexible set of minority opinion holders can 
win over the rest of the population. 

The effect of having un-influencable agents has been 
considered to some extent in prior studies. Biswas et al. 
considered for two state opinion dynamics models 
in one dimension, the case where some individuals are 
“rigid” in both segments of the population, and studied 
the time evolution of the magnetization and the fraction 
of domain walls in the system. Mobilia et. al. (25) con- 
sidered the case of the voter model with some fraction 
of spins representing “zealots” who never change their 
state, and studied the magnetization distribution of the 
system on the complete graph, and in one and two di- 
mensions. Our study differs from these not only due to 
the particular model of opinion dynamics, but also in its 
explicit consideration of different network topologies and 
of finite size networks, specifically deriving how consen- 
sus times scale with network size for the particular case 
of the complete graph. Furthermore, the above men- 
tioned studies do not explicitly consider the initial state 
that we care about - one where the entire minority set 
is un-influencable. A notable exception to the latter is 
the study by Galam and Jacobs [26] in which the authors 
considered the case of “inflexibles” in a two state model of 
opinion dynamics with opinion updates obeying a major- 
ity rule. While this study provides several useful insights 
and is certainly the seminal quantitive attempt at under- 
standing the effect of committed minorities, its analysis 
is restricted to the mean-field case, and has no explicit 
consideration of consensus times for finite systems. 


II. COMPLETE GRAPHS 


A. Infinite network size limit 


We start along similar lines as by considering the 
case where the social network connecting agents is a com- 
plete graph with the size of the network N — oo. We 
designate the densities of uncommitted nodes in states 
A, B as n,,ngp respectively. Consequently, the density of 
nodes in the mixed state AB is nag =1—p—n,— np, 
where p is the fraction of the total number of nodes, that 
are committed. Neglecting correlations between nodes, 
and fluctuations, one can write the following rate equa- 
tions for the evolution of densities: 


dna 2 3 
— = -nNanept+nap+NABNA + TPNAB 
dt 2 
dnp 
= —nang tn yp +Napne — png (1) 


dt 


Before interaction| After interaction 


TABLE I: Shown here are the possible interactions in the bi- 
nary agreement model. Nodes can possess opinion A, B or 
AB, and opinion updates occur through repeated selection 
of speaker-listener pairs. Shown in the left column are the 
opinions of the speaker (first)and listener (second) before the 
interaction, and the opinion voiced by the speaker during the 
interaction is shown above the arrow. The column on right 
shows the states of the speaker-listener pair after the interac- 
tion. 


The terms in these equations are obtained by consider- 
ing all interactions which increase (decrease) the density 
of agents in a particular state and computing the prob- 
ability of that interaction occurring. Table 1 lists all 
possible interactions. As an example, the probability of 
the interaction listed in row eight is equal to the proba- 
bility that a node in state AB is chosen as speaker and a 
node in state B is chosen as listener (nagna) times the 
probability that the speaker voices opinion A (3). 


The fixed-point and stability analyses (see Appendix) 
of these mean-field equations show that for any value of p, 
the consensus state in the committed opinion (n4 = 1—p, 
np = 0) is a stable fixed point of the mean-field dynam- 

2 
ics. However, below p = pe = 3 — 3 (X 54/24 = 1) = 
2 
3 ( V5 —V/24- 1) , two additional fixed points appear: 
one of these is an unstable fixed point (saddle point), 
whereas the second is stable and represents an active 
steady state where n4,ng and nap are all non-zero (ex- 
cept in the trivial case where p = 0) . Fig. [I{a) shows 
(asterisks) the steady state density of nodes in state B 
obtained by numerically integrating the mean-field equa- 
tions at different values of the committed fraction p and 
with initial condition na = 0,ng = 1-p. As p is in- 
creased, the stable density of B nodes ng abruptly jumps 
from % 0.6504 to zero at the critical committed fraction 
pc. A similar abrupt jump also occurs for the stable den- 
sity of A nodes from a value very close to zero below 
Pe, to a value of 1, indicating consensus in the A state 
(not shown). In the study of phase transitions, an “order 
parameter” is a suitable quantity changing (either contin- 


uously or discontinuously) from zero to a non zero-value 
at the critical point. Following this convention, we use 
npg - the density of uncommitted nodes in state B - as 
the order parameter appropriate for our case, character- 
izing the transition from an active steady state to the 
absorbing consensus state. 

In practice, for a complete graph of any finite size, 
consensus is always reached. However, we can still probe 
how the system evolves, conditioned on the system not 
having reached consensus. Fig. [i{a) shows the results 
of simulating the binary agreement model on a complete 
graph for different system sizes (solid lines). For p < pe, 
in each realization of agreement dynamics, neglecting the 
initial transient, the density of nodes in state B, npg, fluc- 
tuates around a non-zero steady state value, until a large 
fluctuation causes the system to escape from this active 
steady state to the consensus state. Fig.[I{a) shows these 
steady state values of ng conditioned on survival, for sev- 
eral values of p. As expected, agreement of simulation 
results with the mean-field curve improves with increas- 
ing system size, since Eq. represents the true evolu- 
tion of the system in the asymptotic network size limit. 
Accordingly the critical value of the committed fraction 
obtained from the mean-field equations is designated as 
p-(oo), although, for brevity, we refer to it simply as pe 
throughout this paper. 

The existence of the transition as p is varied and when 
the initial condition for densities is (n4 = 0,ng = 1 — p) 
can be further understood by observing the motion of the 
fixed points in phase space. Fig. [I{b) shows how the sta- 
ble fixed point and the unstable fixed point move in phase 
space as p is varied from 0 to pe. The active steady state 
moves downward and right while the saddle point moves 
upwards and left. At the critical value pe the two meet 
and the only remaining stable fixed point is the consen- 
sus fixed point. A similar observation was made in the 
model studied in [26]. The fact that the value of ng con- 
verges to © 0.65 and does not smoothly approach zero 
as the stable fixed point and the saddle point approach 
each other, explains the origin of the first-order nature 
of the phase transition. Fig ]shows the representative 
trajectories obtained by integrating the mean-field equa- 
tions for the cases where p = 0.05 (< pe) and p = 0.1 
(> Pe). 


B. Finite network size: Scaling results for 
consensus times 


Even though consensus is always reached for finite N, 
limits on computation time prohibit the investigation of 
the consensus time, Te, for values of p below or very close 
to pe. We therefore adopt a semi-analytical approach pre- 
scribed in that allows us to estimate the consensus 
times for different N for an appreciable range of p includ- 
ing values below pe. We start with the master equation 
which describes the evolution of the probability that the 
network of size N has n (m) uncommited nodes in state 


A (B). We denote by c, the number of committed nodes, 
and by l(= N — n — m — c), the number of uncommitted 
nodes in state AB. 


dpnm 1 1 
dt N N2? 


( — Pnm (2in + Sle +2nm + 1(l— 1) 


1+1 -—1 
oh 2lm + mc) + Pmt UE Fo) 


(n+1)(2m +1-1) 
T Priim 
(1+2)(1+1) 

2 
3(1+.1)(m—1) 


+ Pn-2,m 


T Pnym-1 


(m+ 1)(2n+ 2c+1—-1) 
t Panman 


(1+ 2)(1+1) 
a 


(2) 


T Pn,m—-2 


The factor of 1/N in the LHS comes from the fact that 
a transition between states takes place in an interval of 
time 1/N. The transition rates in each term are the prod- 
uct of two densities, which is responsible for the overall 
factor of 1/N? in the RHS. The probabilities are defined 
over all allowed states of the system (ie. 0 < n < N—c, 
and 0 < m < N—c-—vn for given n ) and the al- 
lowed transitions from any point {n,m} in the interior 
of this state-space are {n,m} > {n,m + 1}, {n,m} > 
{n£1,m}, {n,m} > {n,m + 2}, {n,m} > {n + 2,m}. 


We know from the mean field equations that in the 
asymptotic limit, and below a critical fraction of com- 
mitted agents, there exists a stable fixed point. For fi- 
nite stochastic systems, escape from this fixed point is 
always possible, and therefore it is termed metastable. 
For a finite system, the probability of having escaped 
to the metastable fixed point as a function of time is 
P.(t) = 1 — P,(t) where P,(t) is the survival probabil- 
ity. The surviving fraction is constrained to be in the 
allowed region of n,m quadrant excluding the true fixed 
point {N —c,0}. If the number of committed agents is 
far lower than pe N we expect that this surviving frac- 
tion will occupy configurations around the metastable 
fixed point, and the occupation probabilities Pn,m will 
be peaked around the metastable fixed point. In sys- 
tems which exhibit such long lived metastable states in 
addition to an absorbing fixed point, applying a quasista- 
tionary (QS) approximation has been found to be useful 
in computing quantities of interest RH29). This approx- 
imation assumes that after a short transient, the occu- 
pation probability, conditioned on survival, of allowed 
states excluding the consensus state, is stationary. Fol- 
lowing this approximation, the distribution of occupation 
probabilities conditioned on survival can be written as, 
Düm. = Pnm(t)/Ps(t) and using this form in the mas- 
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FIG. 1: (a) The steady state density ng of nodes in state B as a function of committed fraction p for complete graphs of 
different sizes, conditioned on survival of the system. Simulation results are from 100 realizations of the binary agreement 
dynamics. (b) Movement of the stable fixed point and the saddle point in phase space as a function of committed fraction p 
(see text). The location of these points in phase space was obtained through a stability analysis (not shown) of the mean field 
equations Eq. (I). The inset shows the density of nodes in state B at the stable (red) and unstable(blue) fixed points as p is 
varied; the expressions from which these numerical values are obtained are given in the Appendix. 


ter equation (Eq. (2)), we get: 


dP,(t) ~ 


P(t) 
dt Pram = 


=n (Ban (2ln + se + 2nm 


+ I(l- 1)+2lm + me) 
X 3(1+1)(n—1+c) 
Pn-im ~p 


2 (n+ 1)(2m+1—- 1) 
T Prim 


4 (1+ 2)(1+1) 
EE Pn-2m z 


3(1+1)(m-1) 


— Pnym-1 


> (m +1)(2n +2c+1l-— 1) 
~ a, 


Z Cti) (3) 


— Pn,m—2 2 


Considering transitions from states {N — c — 1,0} and 
{N —c— 2,0} to the absorbing state {N — c, 0}, we obtain 
the decay rate of the survival probability dP,(t)/dt: 


oo = —P,(t) för- Hi + Pr-o-20 


(4) 
Substituting Eq. (4) into Eq. (8), we finally obtain a con- 
dition that the occupation probabilities conditioned on 
survival must satisfy : 


5 Qnm 
nm = o = 5 


where, nm = Qnm(t)/Ps(t) is obtained through explicit 


consideration of the terms in the master equation: 


~ R 3(1+1)(n-—l+e 
Qnm = a EU 


2 
5 n+1)\(Q2m+1-1 7 1+ 2)? 
Tr Freee a X ) fy ee a 
2 2 
` 3(1+1)(m-1 B t+2y 
T foa UE oS ) 
2 2 
z m+ 1)(2n+2c+l-—1 
T Ptit M 2 ) (6) 
and Qo = [Pn-e-1,0 (= +2 By-c-20| is the 


term arising from the decay of the survival probability 
(Eq. (2). Wnm is the coefficient of pram (Pnm) within 
the brackets on the right hand side of Eq. (2) (Eq. BJ) 
and is equal to the transition rate out of state {n,m} 
times N?. 

Eq. indicates that the survival probability decays 
exponentially with a rate A = Qo /N. Since the mean life- 
time of an exponentially decaying process is the inverse 
of the decay rate, it follows that the mean consensus time 
(neglecting the short transient before the QS distribution 
is attained) is 


1 i aNEDY .. 2 
Te x—-=1 —c— Ar? —c—2,057 
X / [øx 1,0 ( aN ) + PN som 


(7) 

Thus, knowledge of nms (in particular, py_—-—1,9 and 
PN—c—2,0) would allow us to calculate Te through Eq. (7). 
In order to obtain ppm (for all 0 <_n,m), we adopt 
the iterative procedure proposed in . Following this 
procedure, we start with an arbitrary initial distribu- 


tion p? m, and obtain a new distribution using: pit) = 


a Pom +(1—a) Pam 


Wi Op where 0 < a < 1 is an arbitrary 
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FIG. 2: Trajectories (obtained from integration of the mean-field equations Eq. (I)) in the phase-plane show the nature of flows 
from different regions of the phase-plane into existing fixed points for a) p = 0.05 (< p-) and b) p=0.1 (> pc) 


parameter, and Ĝim, Wf, and Qj, are all obtained using 
the probability distribution at the current iteration, Pim 

With sufficient number of iterations, this procedure 
is expected to converge to a distribution that satisfies 
Eq. and which is thus, the desired QS distribution. 
In our case, we obtained acceptable convergence with a 
choice of a = 0.5 and 30000 iterations. 


Following the above method, we obtain the QS distri- 
bution, and consequently the mean consensus times Te 
for different values of committed fraction p and system 
size N. Fig. B{a) shows how the consensus time grows 
as p is decreased beyond the asymptotic critical point pe 
for finite N. For p < pe, the growth of Te is exponential 
in N (Fig. Bib), consistent with what is known regarding 
escape times from metastable states. For p > pe, the QS 
approximation does not reliably provide information on 
mean consensus times, since consensus times themselves 
are small and comparable to transient times required to 
establish a QS state. However, simulation results show 
that above pe the scaling of mean consensus time with 
N is logarithmic (Fig. B{c)). A snapshot of the QS dis- 
tribution (Fig. Æ) near pe (p = 0.09) for a system of size 
N = 100 shows clearly the bimodal nature of the distri- 
bution, with the two modes centered around the stable 
fixed point, and the consensus fixed point. 


The precise dependence of consensus times on p can 
also be obtained for p < pe by considering the rate of ex- 
ponential growth of Te with N. In other words, assuming 
T. ~ exp(a(p)N), we can obtain a(p) as a function of p. 
Figure Bd) shows that a(p) ~ |p — pe|” where v ~ 1.65. 
Thus, below pe, we have: 


Te(p < pe) ~ exp((pe — p)” N) (8) 


This exponential growth is presumably modulated by fac- 
tors of log N which become dominant only when p = pe. 
Above pe, the dependance of Te on p as seen from simu- 
lations is negligible (not shown). 


II. SPARSE NETWORKS 


Next, we present simulation results for the case when 
the underlying network topology is chosen from an en- 
semble of Erdős-Rényi (ER) random graphs with given 
size N and given average degree (k). The qualitative fea- 
tures of the evolution of the system in this case are the 
same as that of the complete graph, although the crit- 
ical fraction pe displays some dependence on (k). For 
small (k) and fixed N, the drop in consensus times oc- 
curs slightly earlier in p for ER graphs than for a complete 
graph of the same size, as shown in Fig. [b{a). However 
for p > pe, a complete graph has shorter consensus times 
(on average) than an ER graph of the same size. Above 
Pe, the difference between consensus times for a graph 
with average degree (k) and the complete graph, AT,, 
decays approximately as power law with increasing (k) 
(Fig. BYb)). The deviation from a perfect power law is 
likely due to other weaker (k} dependent terms, presum- 
ably logarithmic in (k). We also performed simulations 
of the binary agreement model on Barabasi-Albert (BA) 
networks (Fig. BẸ)), and found similar qualitative behav- 
ior as observed for ER networks including the difference 
from mean-field behavior. We leave a detailed analysis of 
the dependence of the critical fraction pe and the consen- 
sus times T, on the average degree (k) of sparse networks, 
for future work. 


IV. SUMMARY 


In closing, we have demonstrated here the existence of 
a tipping point at which the initial majority opinion of a 
network switches quickly to that of a consistent and in- 
flexible minority. There are several historical precedents 
for such events, for example, the suffragette movement 
in the early 20th century and the rise of the American 
civil-rights movement that started shortly after the size of 
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FIG. 3: (a) Mean consensus time Te for p < pe obtained by using the QS approximation. (b) Exponential scaling of Te with N, 
for p < pc; mean consensus times (circles, squares) are obtained using the QS approximation. The lines are guides to the eye. 
(c) Logarithmic scaling of Te with N for p = 0.3 > pe; mean consensus times are obtained from simulations. The line shows 
the best linear fit to the data. (d) The rate, a(p), of exponential growth of the consensus time with N as a function of p — pe 
(see text). Circles show the values of a(p) obtained for p = 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8 by considering the scaling of Te with 
N, for these values of p. The straight line shows a linear fit to the data plotted on a log-log scale. 


the African-American population crossed the 10% mark. 
Such processes have received some attention in sociologi- 
cal literature under the term minority influence bod Bi. 
Our motivation here has been to study this process in 
more detail through semi-analytical methods and simu- 
lations for finite-sized and sparse networks, within the 
realm of a particular social influence model - the binary 
agreement model. There are several open questions and 
extensions of this work that are worth studying, in our 
opinion: for example, given a network with non-trivial 
community structure, what is the optimal scheme for se- 
lecting committed agents (for a given committed frac- 
tion) that would minimize consensus times, and reduce 
Pe? Secondly, extensions of the model to include utility- 
driven opinion switching by agents may be useful in de- 
signing optimal incentive schemes for opinion spreading. 
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FIG. 4: The quasistationary distribution nm for p = 0.09 
and N = 100. 


Appendix: Fixed points of the mean-field equations 


Here, we analyze the mean-field equations for the ex- 
istence of fixed points. To simplify notation we use the 
notation x = na and y = ng. Thus for a fixed point of 
the evolution given by Eq. (I): 


-zy + (l-g-y-p)?+2(1-x2—y-p) 
+ 1.5p(1-—a2-—y-p)=0 
-zy + (L-x-y-p)? +y(l-zr-y-p) 
— yp=0 (9) 


which can be reduced further to: 


a = [(1—y—p/4)* — 9p?/16]/(p/2 + 1) 
y = (l-w—p)? (10) 


Substituting the expression of x into the expression for 
y, and denoting z? = y we get 


2(z? — (2—p/2)z-+p/2+1)=0 (11) 


For any value of p, z = zp = 0 is a solution to the above 
equation. In other words, for any value of p, the mean- 
field equations admit a stable fixed point, ng = z1 = 
1 — p, ng = yı = 0 which represents the network having 
reached a consensus state where all nodes have adopted 
the opinion of the committed agents. 

The remaining fixed points are roots of 


f(z) = 2 —(2—p/2)z+-p/2+1=0. (12) 


In order to find the criterion which has to be satisfied 
for valid roots (ie. 0 < z < \/((1 —>p))) of the above 
equation to exist, we analyze the extrema of the function 
f(z) which are given by: 


f'(z) = 32? —2+7/2=0 (13) 
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FIG. 5: (a) The steady state density ng of nodes in state B 
as a function of committed fraction p for Erdés-Rényi graphs 
of different sizes with (k) = 10, conditioned on survival of the 
system. Lines are mean values of ng obtained from simula- 
tions; asterisks show mean-field consensus times obtained by 
integrating Eq. (b) Scaling of AT.((k)) (defined in text) 
with (k); the line shows the best linear fit to the data. (c) 
The steady state density ng of nodes in state B as a function 
of committed fraction p for Barabasi-Albert (BA) networks 
of different sizes with (k) = 10, conditioned on survival of the 
system 


Thus, the extrema occur at : 


21,2 = +y 2/3 — p/6 (14) 


It can be seen from Eq. that f(z) is increas- 
ing, decreasing and increasing again in the intervals 
(—00, z21),(21, 22),(22, +00) respectively. Consequently, 


f(z) achieves a maximum at —1 < z1 = —,/2/3—p/6 < 
0 and a minimum at 0 < z2 = \/2/3 — p/6 < 1. Further- 
more, since f(—2) = —p/2—3 < 0 and f(—1) = 2 > 0, 
one root of f(z) = 0 occurs in the interval —2 < z < —1. 
Since f(z) is positive at z1, decreasing from 21 to z2 where 
it achieves a minimum, and increasing thereafter, it fol- 
lows that a necessary and sufficient condition for more 
roots of f(z) = 0 to exist, is that f (z2) be less than zero: 


f (2) = 23 — (2 — p/2)z2 + p/2+1<0 


Denoting z2 = q and p = 4 — 6q? (from Eq. (14)) yields 
the following inequality for q as a condition for more roots 
of f(z) = 0 to exist: 


fla) =@P41.5¢-15>0 (15) 


(Note that z2 is itself a function of p). Analyzing the 
derivative of f(q) enables us to glean that the inequality 
Eq. is satisfied for q > qo where qo is the solution of 
the cubic equation f(q) = 0 and is given by: 


qo = [i/s4 vai s-vai-1| /2 


Thus, the original fixed point equation Eq. (9) has at 
least one valid root besides z = 0, so long as p is less 
than or equal to: 


5 3/s e Bhs 
Peo==—=( 5424-1 =5 5—V/24-1 


2 2 

(16) 
which using standard computer algebra software is eval- 
uated to be pe = 0.09789. Using, z? = y = qo and 
Eq. (£0), we obtain the state of the system at pe to be 
{na, np} = {0.0957, 0.6504}. It also follows from the ex- 
pression for f(z), that f(0) > 0 and therefore if f(z2) is 
negative, Eq. has two roots on the positive line when 
p < pc. Thus there are two fixed points of Eq. (9) when 
P < Pe. 

The exact expressions for these fixed points (that can 
also be obtained numerically), obscure their dependence 
on p. We therefore adopt an approximation which ex- 
hibits a much clearer dependence of the fixed point values 
on p, and numerically yield values close to those obtained 
from the exact expressions. Substituting z = t,/2 — a in 
Eq. reduces it to: 


2 


B—t+r=0 (17) 
where r = ace Clearly, r is a monotonically 


increasing function of p, and therefore nA <re< 


ltpe/2 0 _ 2 Sela 
(/2—pef22 «BV for 0 < p < pe, our range of inter 


est. Function g(t) = t? — t is monotonically decreasing 


for t < —1 and g(—1) = 0, while g(—2) < —1 < —-r. 
Hence, there is a real root tı € (—2,—1) to Eq. (7) 
which is a monotonically decreasing function of r, but 
which clearly does not yield a valid fixed point. This 
root can be expressed as t1(r) = = + e(r), where e(r) 
is monotonically decreasing from less than 0.0106 to 0 
over the range of our interest for r. Substituting this ex- 
pression back into Eq. and neglecting powers of e(r) 
higher than unity, we get an approximation of € in terms 
of r, and consequently an approximation for tı: 


16 
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with relative error of less than 0.01%. Now, we can fac- 
torize the LHS of Eq. (17%) and write it as (t? +bt+c)(t— 
tı). Equating this factorized expression with t — t +r, 
gives us b and c in terms of r. Thus, two more roots 
of Eq. @7) are obtained in terms of r by solving the 
quadratic equation t? + bt + c = 0 which yields: 


ti(r) & (18) 


wir 


8 r 
t23 = —= + = E4] >- sez SS 19 
TT 81 (13) 


Finally, we can obtain the values of z associated with the 
above roots, and therefore the values of x and y written 
in terms of these roots are derived as: 
4-p 
ns = Ba E 
(4 — 4y2,3 — p)? — 9p? 


— 4y- p) — Ip 2 
v23 8p + 16 (20) 


The stability of these fixed points can be checked via 
linear stability analysis. Following the standard proce- 
dure, the stability matrix is given by: 


g- -1- $ —2 + 2y* +5 
—2 + 2x* + 2p —1 


where (x*, y*) is the fixed point under consideration. The 
eigenvalues of the stability matrix are given by: 


1 


à= i(-4-p+ (17p? + 64(2* — 1)(y* — 1) 


+ 16p(a* + 4y* — st Jan 


From the expression for the eigenvalues we numerically 
determine that the real part of both the eigenvalues is 
negative for (x2, y2) over the range 0 < p < pe indicating 
that (x2, y2) that is a stable fixed point. This is however, 
not the case for (x3, y3), making it unstable. Similarly, 
the consensus fixed point (xı, yı) is found to be stable 
for 0 < p < 1. Finally, we note that as p —> 0, the 
stable fixed point converges to na = 0,ng = 1, while the 
unstable fixed point converges to na = npg & 0.38. 
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